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1 Introduction 

Understanding the statics and dynamics of elastic systems in a random en- 
vironment is a longstanding problem with important applications for a host 
of experimental systems. Such problems can be split into two broad cate- 
gories: (i) propagating interfaces such as magnetic [1, 2, 3, 4] or ferroelectric 
[5, 6] domain walls, fluid invasion in porous media [7], contact line in wetting 
[8], epitaxial growth [9] or crack propagation [10]; (ii) periodic systems such 
as vortex lattices [11, 12, 13], charge density waves [14], or Wigner crystals 
of electrons [15, 16]. In all these systems the basic physical ingredients are 
identical: the elastic forces tend to keep the structure ordered (flat for an 
interface and periodic for lattices), whereas the impurities locally promote 
the wandering. From the competition between disorder and elasticity emerges 
a complicated energy landscape with many metastable states. This results 
in glassy properties such as hysteresis and history dependence of the static 
configuration. 

To study the statics, the standard tools of statistical mechanics could be 
applied, leading to a good understanding of the physical properties. Scaling 
arguments and simplified models showed that even in the limit of weak dis- 
order, the equilibrium large scale properties of disordered elastic systems are 
governed by the presence of impurities. In particular, below four (internal) 
dimensions, displacements grow unboundedly [17] with the distance, resulting 
in rough interfaces and loss of strict translational order in periodic structures. 
To go beyond simple scaling arguments and obtain a more detailed description 
of the system is rather difficult and requires sophisticated approaches such as 
replica theory [18] or functional renormalization group [19]. Much progress 
was recently accomplished both due to analytical and numerical advances. 
For interfaces, the glassy nature of the system is confirmed (so called random 
manifold), and a coherent picture of the system is derived from the various 
methods. Periodic systems have also been shown to have glassy properties but 
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to belong to a different universality class than interfaces, with quite different 
behavior for the long range nature of the correlation functions [20, 12, 13]. 

The competition between disorder and elasticity manifests also in the dy- 
namics of such systems, and if any in a more dramatic manner. Among the 
dynamical properties, the response of the system to an external force F is 
specially crucial, both from a theoretical point of view, but also in connec- 
tion with measurements. Indeed in most systems the velocity v versus force 
F characteristics is directly measurable and is simply linked to the transport 
properties (voltage-current for vortices, current-voltage for CDW and Wigner 
crystals, velocity-applied magnetic field for magnetic domain walls). 

Some of the questions related to this issue are shown in Fig. 1. In the pres- 




Fig. 1. The velocity v induced by an external force F of a disordered elastic 
system. At zero temperature T — the system stays pinned until a critical force F c 
is reached. At finite temperature a motion can occur even for forces below threshold 
F < F c since the barriers to motion can always be passed by thermal activation. 



ence of disorder it is natural to expect that, at zero temperature, the system 
remains pinned and only polarizes under the action of a small applied force, 
i.e. moves until it locks on a local minimum of the tilted energy landscape. At 
larger drive, the system follows the force F and acquires a non-zero asymp- 
totic velocity v. So a first set of questions is prompted by the zero temperature 
properties. What is F c ? An estimate of F c can be obtained via scaling argu- 
ments [21] or with a criterion for the breakdown of the large velocity expansion 
[22, 23] and related to static quantities such as the Larkin-Ovchinikov length, 
or computed numerically by an exact algorithm [24] . The v — F curve at T = 
is reminiscent of the one of an order parameter in a second order phase transi- 
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tion [25] . Here the system is out of equilibrium so no direct analogy is possible 
but this suggests that one could expect v ~ (F — F C )P with a dynamical crit- 
ical exponent (3 [26]. Calculation of such exponents is of course an important 
question [27, 28, 29, 30]. 

Another important set of questions pertain to the nature of the moving 
phase itself, and in particular to the behavior at large velocity: how much this 
moving system resembles or not the static one [31, 32, 33] ? This concerns 
both the positional order properties and the fluctuations in velocity such as 
the ones measured in noise experiments. 

Finally, one of the most important questions, and the one on which we 
will concentrate in these notes, is the response well below threshold F <C F c 
at finite temperature. In this regime, the system is expected to move through 
thermal activation. What is the nature of this motion and what is the velocity? 
The simplest answer would be that the system can overcome barriers via 
thermal activation, [34] leading to a linear response at small force of the form 
- U >/ T F, where U b is some typical barrier. However it was realized that 
such a typical barrier does not exist in a glassy system [35, 36, 37, 38] and 
that the response of the system was more complicated. The motion is actually 
dominated by barriers which diverge as the drive F goes to zero, and thus 
the flow formula with finite barriers is incorrect. Well below F c , the barriers 
are very high and thus the motion, usually called "creep" is extremely slow. 
Scaling arguments, relying on strong assumptions such as the scaling of energy 
barriers and the use of statics properties to describe an out of equilibrium 
system, were used to infer the small F response. This led to a non linear 
response, characteristic of the creep regime, of the form v ~ exp(— C F~^/T). 

Given the phcnomcnological aspect of these predictions and the uncon- 
trolled nature of the assumptions made, many open questions remain to be 
answered, in particular whether such a behavior is indeed correct [39] and can 
be derived directly from microscopic equation of the motion [40, 29]. We will 
review these questions in these notes. The plan of the notes is as follows. In 
Sec. 2 we recall the basic concepts of interfaces in the presence of disorder. In 
Sec. 3, we recall the phenomenological derivation of the creep law. We present 
the microscopic derivation of the creep law from the equations of motion in 
Sec. 4, and discuss the similarities and differences with the phenomenological 
result. In Sec. 5 we focus on the situation of domain walls. Such a situation is 
a particularly important both for experimental realizations of the creep but 
also because one dimension is the extreme case for such systems. Conclusions 
can be found in Sec. 6. 



2 Basic concepts 

Let us introduce in this section the basic ingredients of the systems under 
study. We will focus in these short notes to the case of interfaces, but similar 



4 T. Giamarchi, A. B. Kolton, and A. Rosso 



concepts apply to periodic systems as well. The interface is a sheet of dimen- 
sion d living in a space of dimensions D. For realistic interfaces D = d + 1 
but generalization are of course possible (for example D = d corresponds to 
periodic systems). We call r the internal coordinate of the interface and z all 
its transverse directions. The interface position is labelled by a displacement 
it(r) from a flat configuration. This determines totally the shape of the inter- 
face provided that u is univalued, i.e. that there are no overhangs or bubbles. 
The case of a one dimensional interface (d = 1) in a two dimensional film 
is shown in Fig. 2. Since the interface distortions cost elastic energy, its zero 



Fig. 2. A one dimensional interface (domain wall) living in a two dimensional space 
(film) . The position of the interface is determined (provided there are no overhangs 
or bubbles) by the displacement u from a flat configuration. 

temperature equilibrium configuration in the absence of disorder is the flat 
one. Deviation from this equilibrium position are described by an Hamiltonian 
H [it] which is a function of the displacements u. For small displacements one 
can make the usual elastic approximation 



where u q is the Fourier transform of u(r) and c(q) are the so called elastic 
coefficients. If the elastic forces acting on the interface are short ranged then 
one has c(q) = cq 2 which corresponds to 
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For some interfaces where long range interactions play a role different forms 
for the elasticity are possible. This is in particular the case when dipolar forces 
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[41] are taken into account [6] or for the contact line in wetting [42] and crack 
propagation [43]. 

In addition to the elastic energy the interface gains some energy by cou- 
pling to the disorder. Two universality classes for the disorder exist (see 
Fig. 3). The so called random bond disorder corresponds to impurities that 
directly attract or repel the interface. On the contrary, for the so called ran- 
dom field disorder the pinning energy is affected by all the randomness that 
the interface has encountered in its previous motion. If V(z,r) denotes the 
random potential generated by the impurities the pinning energy writes: 
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Fig. 3. The two types of disorder (the names are coming from the magnetic real- 
ization of such systems). The dark circles are the impurities that contribute to the 
pinning energy of the interface. In the random bond case only neighboring impuri- 
ties contribute while in the random filed case all the impurities on the left side of 
the the interface contribute. This makes the latter disorder effectively long ranged, 
even if the disorder potential V(z, r) is short range. 



The competition between disorder and elasticity manifests itself in the 
static properties of the interface. The presence of disorder leads to the ap- 
pearance of many metastable states and glassy properties. In particular, the 
interface deviates from the flat configuration and becomes rough. From the 
scaling of the relative displacements correlation function, a roughness expo- 
nent £ can be defined by 

B(r) = (Kr) -u(0)] 2 } oc r 2C (4) 

where ( ) denotes thermodynamic average and 7TT denotes disorder average. 
We will not enter in more details about the statics here and refer the reader 
to the literature on that point [44, 9]. 
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Dynamics is much more complicated since the standard tools of statistical 
physics can not be used. One has to study the equation of motion of the 
system 

This equation is written for overdamped dynamics, but can include inertia 
as well. r\ is the friction taking into account the dissipation, F the external 
applied force, and C( r ?0 a thermal noise, needed to reproduce the effect of 
finite temperature. The correlation of the thermal noise is (£(r, t)C,(r' , t')) — 
2rjT5(r — r')S(t — t 1 ). Solving this equation of motion allows to extract all the 
dynamical properties of the system. The presence of disorder in the Hamilto- 
nian H makes this a very complicated proposal. In the absence of the external 
force F = 0, this Langevin equation allows to recover the static properties af- 
ter the system has achieved its thermal equilibrium. 



3 Creep, phenomenology 

Let us focus here on the response of the system to a very small external 
force. For usual systems we expect the response to be linear. Indeed earlier 
theories of such a motion found a linear response. The idea is to consider 
that a blob of pinned material has to move in an energy landscape with 
characteristic barriers Ub as shown in Fig. 4. The external force F tilts the 



energy landscape 




Fig. 4. In the Thermally Assisted Flux Flow (TAFF) [34] a region of pinned 
material is considered as a particle moving in an energy landscape characterized by 
characteristic barriers Ub- This leads to an exponentionally small but linear response. 



energy landscape making forward motion possible. The barriers are overcomed 
by thermal activation (hence the name: Thermally Assisted Flux Flow) with 
an Arrhenius law. If the minima are separated by a distance a the velocity is 

V OC e-m-F*M - e -P(U b +Fa/2) „ e -0U b p (g) 

The response is thus linear, but exponentially small. 



Dynamics of disordered elastic systems 



7 



However this argument is grossly inadequate for a glassy system. The rea- 
son is easy to understand if one remembers that the static system is in a 
vitreous state. In such states a characteristic barrier U h does not exist, since 
barriers are expected to diverge as one gets closer to the ground state of the 
system. The TAFF formula is thus valid in systems where the glassy aspect is 
somehow killed and the barriers do saturate. This could be the case for exam- 
ple for a finite size interface. When the glassy nature of the system persists up 
to arbitrarily large length scales the theory should be accommodated to take 
into account the divergent barriers. This can be done quantitatively within 
the framework of the elastic description [35, 36, 38, 45]. The basic idea rests 
on two quite strong but reasonable assumptions: (i) the motion is so slow that 
one can consider at each stage the interface as motionless and use its static 
description; (ii) the scaling for barriers, which is quite difficult to determine, 
is the same as the scaling of the minimum of energy (metastable states) that 
can be extracted again from the static calculation. If the displacements scale 
as u ~ L'' then the energy of the metastable states (see (2)) scales as 

E(L) ~ L d - 2+2 < (7) 

where we use that elastic and pinning energy scale the same way. Since the 
motion is very slow, the effect of the external force is just to tilt the energy 
landscape 

E(L) -F J d d ru(r) ~ L d - 2+2<: - FL d+ ^ (8) 

Thus, in order to make the motion to the next metastable state, one needs to 
move a piece of the pinned system of size 

Lopt-Q;) 2 " 4 (9) 

The size of the optimal nucleus able to move thus grows as the force decrease. 
Since the barriers to overcome grow with the size of the object, the minimum 
barrier to overcome {assuming that the scaling of the barriers is also given by 
(7)) 

d-2+2C 

" (10) 
leading to the well known creep formula for the velocity 



v cx exp 



-We ( f 



F \ M 1 

(11) 



where F c is the depinning force and U c a characteristic energy scale and the 
creep exponent \i is given by, 
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Equations (11) and (12) are quite remarkable. They relate a dynamical prop- 
erty to static exponents, and shows clearly the glassy nature of the system. 
The corresponding motion has been called creep since it is a sub-linear re- 
sponse. Of course the derivation given here is phenomenological, and it will 
be important to check by other means whether the results here hold. This will 
be the goals of the two next sections, where first the creep law will be derived 
directly from the equation of motion in d = 4 — e dimensions, and then the 
creep will be examined in the important case of d = 1 domain walls. 



4 Around four dimensions 



The previous phenomenological derivation of the creep formula rests on very 
strong hypothesis. In particular it is assumed that: (a) the motion is domi- 
nated by the typical barriers, and not by tails of distributions in the waiting 
times or barriers; (b) the motion is so slow that the line has the time to 
completely re-equilibrate between two hopping events so that one can take 
all exponents as the equilibrium ones. Given the phenomenological ground of 
these predictions and the uncontrolled nature of the assumptions made, both 
for the creep and for the depinning, it is important to derive this behavior in 
a systematic way, directly from the equation of motion. 

In principle one has simply to solve the equation of motion (5). In prac- 
tice this is of course quite complicated. A natural framework for computing 
perturbation theory in off-equilibrium systems is the dynamical formalism 
[46, 47]. Integrating on all configurations u we can exponentiate the equation 
of motion by introducing an auxiliary field u : 



SH 



F-((r,t) 



5u(r, t) 

j VuVu exp 



, du SH 



(13) 



the thermal and disorder average can easily be done, leading to a field theory 
with some action S 



lUrtiilrt ~ F I lU r t 
rt J rt 



S(u, u) = iurtivdt - cV 2 )u rt -rjT 

Jrt Jrt 

— - I iu rt iu rt , A(u rt - u rt ,) (14) 

2 Jrtt' 

where A is defined in the correlator of the pinning force F p = —SH^is/Su, as 
F p (u, r)F p (u', r') = A(u - u)S(r - r'). (15) 



The functional form of this correlator depends on whether one has random 
bond or random field disorder (see e.g. [29] for more details). Essentially A is a 
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function with a width of the order of the correlation of the disorder along the z 
direction. The advantages of this representation are many: disorder and ther- 
mal averages (A[u]) = (.A [it]) 5 of any observable A[u] can be computed with 
the weight e~ s ; the response functions to an external perturbation h rt are 
simply given by correlations with the response field: {A[u]iu rt } = g£—(A[u]}. 
In addition, since we are much more familiar with fields theories than with 
equations of motion, we have at our disposal a variety of tools to tackle the 
action S. Although it is impossible to solve the action exactly it is possible to 




Fig. 5. Although the correlator of the disorder is initially an analytic function, a 
non analyticity (cusp) appears at a finite scale l c . This length scale corresponds to 
the Larkin-Ovchinikov length at which pinning and metastability occur. 



look at its properties using a renormalization group procedure. We will not de- 
tail the procedure but just recall here the resulting functional renormalization 
group (FRG) flow equations [40, 29] to give a flavor of their physics 



dA(u) = (e - 2()A(u) + (uA'(u) + fA"(u) 

+ fe- s - s ' \A"(u) (A((s' - s)X) - A(u + (s' 
Js>o,s'><:r- * 

-A'(u- s'X)A'(u + sX) 
+A'({s' + s)A) (A'(u - s'X) - A'{u + sX) 



(16) 



s)X) 



<91nA = 2-C 



e- s sA"(sX) 



s>0 



<91nf = e-2-2C+ f e - s sXA'"{sX) 

Js>0 



OF = e-( 2 -WcA 2 



e- s A'{sX) 



s>0 



where e = 4 — d, d denotes and A = rji v. The tilde denotes rescaled quan- 
tities (see [29] for the notations). Contrarily to the standard case of critical 
phenomena, where the potential A can be expanded in powers of the field and 
only the first terms are relevant, here all the powers of the expansion have the 
same dimension. It is thus necessary to renormalize the whole function [48] 
(in other words one has an infinite set of coupled renormalization equations). 
One of the most important consequences is shown in Fig. 5: the renormalized 
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function A becomes nonanalytic beyond a certain length scale, and develop 
a cusp which signals pinning and the glassy properties of the system. This 
cusp appears at a finite length scale corresponding to the Larkin-Ovchinikov 
length [23] and it is directly related to the existence of the finite critical force 
F c [27, 28] at zero temperature. The FRG procedure has been push up to two 
loop expansion [49, 50]. 

The presence of a finite temperature and a finite velocity (proportional 
to A) prevents the appearance of the cusp [40, 29]. For very small external 
force, the way the cusp is cut occurs in two steps, as shown in Fig. 6. At 



thermal depinning flat 




Fig. 6. The cusp is rounded by both the temperature and the finite velocity of 
the interface. In the thermal regime, the main source of rounding comes from the 
temperature and the role of velocity is negligible. Then the system enters a regime 
in which the main source of rounding is the velocity and the role of temperature is 
negligible. This regime is very similar to a depinning regime. Finally the velocity 
rounds the whole correlator of the disorder and thus disorder is washed out by the 
averaging due to motion. At this length scale one recovers a purely thermal interface. 

very low velocity the cusp is cut first by the temperature and the velocity 
can be forgotten. A physical way to interpret this regime is that the motion 
consists essentially in overcoming the barriers by thermal activation. This 
regime corresponds essentially to the one in the phenomenological derivation 
of the creep. Increasing the scale I the temperature renormalizes down and 
the velocity renormalizes up. For this reason, above a lengthscale that can 
be identified with L opt , the cusp starts to be regularized by the velocity. The 
temperature can now be forgotten and a regime very similar to the standard 
depinning regime at T = takes place. Then, finally, at a certain length scale 
L v the whole u dependence of the correlator of the disorder A(u) is erased 
by the finite velocity. This corresponds to a regime where the motion of the 
interface has averaged over the disorder and thus, in the moving frame, the 
interface is now simply submitted to the thermal- like noise [27] . 



Dynamics of disordered elastic systems 11 

The FRG calculation of the velocity confirms the phenomenological argu- 
ments and finds the creep law (11). Moreover, at the first order in e, the creep 
exponent /x agrees with the scaling (12). The velocity behavior is thus dom- 
inated by the first (thermal) regime. On the other hand because the second 
regime exists, we expect that the phenomenological derivation is incorrect as 
far as the characteristic lengthscales of the problem are concerned. Indeed, 
the phenomenological derivation would predict that the characteristic size of 
a moving domain is the optimal length scale L pt, which coincides with the 
end of the thermal regime. However the FRG equations predict that the size 
of a moving domain corresponds to L v , the scale of the end of the depinning 
regime. A physical way to understand this behavior is to say that the mo- 
tion of the thermal nucleus, of size L opt , triggers an avalanche of a larger size 
L v . The FRG thus predicts much larger avalanche scales than what should 
be naively expected from the phenomenological theory, as shown in Fig. 7. 
Such large avalanches are in agreement with recent experiments in magnetic 
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Fig. 7. The phenomenological theory of creep (left) would predict that the 
avalanche size corresponds to the size of the thermal nucleus L opt . From the FRG a 
quite different size emerges corresponding to a larger avalanche of size L v triggered 
by the thermal motion, in a way similar than for a depinning process. A and a are 
two characteristic exponents. 



systems [3]. 

Confirming the stretched exponential behavior of the creep is of course 
an experimental challenge given the large span of velocities needed. The first 
unambiguous determination of the creep law with a precise determination of 
the exponent was made in magnetic films, for one dimensional domain walls 
[1], and confirmed with subsequent measurements [3, 4]. Ferroelectric systems 
[5, 6] have shown a creep exponent compatible with two dimensional domain 
walls in presence of dipolar forces . In periodic systems, such as vortices, it is 
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more delicate to determine the precise value of the exponent, even when non 
linear behavior has been clearly observed. The experiment [51] shows a creep 
exponent in agreement with the theoretical predictions. 




U 



Fig. 8. Discretization scheme for the elastic line driven in a random potential. 



5 Low dimensional situation: domain walls 

Around four dimensions the creep hypothesis gives a velocity dependence 
consistent with the one obtained from the microscopic derivation, at least up 
to the order e at which the renormalization group analysis can be performed. 
Let us now focus at the other extreme limit, namely when the wall is one 
dimensional and moves in a two dimensional space. The interest in such a 
situation is twofold. First, as already mentioned, controlled experiments are 
performed on domain wall motion. Second, from the theoretical point of view 
the situations of a low dimensional domain wall is very interesting. Thermal 
effects are increasingly important as the dimension is lowered. For d < 2 they 
lead to a roughening of the domain wall, even in the absence of disorder (with 
an exponent £t = (2 — <£)/2). One can thus expect more intricate competition 
between temperature and disorder effects. 

Numerical simulations are a valuable alternative theoretical tool to address 
this open issue. In this respect, Langevin dynamics simulations have been 
used to study both the velocity-force (v-F) characteristics and the dynamic 
roughness ( of an elastic string in a random potential [52, 53, 39]. In [39] we 
have studied equation (5) with a short range elasticity: 



Dynamics of disordered elastic systems 



13 



77-u(r, t) = cd 2 r u{r, t) + F p (u, r) + F + C(r, t) (17) 

where F p (u,r) — —d u V(u(r),r) is the pinning force derived from the random 
bond disorder V(u, r). 

To solve numerically (17) we discretize the string along the r direction, 
r — > j = 0, . . . , L — 1, keeping Uj(t) as a continuous variable. A second order 
stochastic Runge-Kutta method [54] is used to integrate the resulting equa- 
tions. To model a continuous random potential we generate, for each j, a 
cubic spline V(uj,j) passing through regularly spaced uncorrelated Gaussian 
random points [24]. The geometry of our system is shown in Fig. 8. We are 
interested in the v-F characteristics. Typical curves, obtained in the simula- 
tions, are shown in Fig. 9. In the whole range of temperature and pinning 




Fig. 9. After [39]. v-F characteristics for several temperatures, increasing from 
bottom to top. Solid lines are fits of the creep formula (11) with U c and n as fitting 
parameters. The inset assures the validity of the creep formula in the range of 
temperature and velocity analyzed. 



strength analyzed we find that the v-F curve can be well fitted by the creep 
formula (11) with U c and /i as fitting parameters. We thus confirm the pre- 
dicted stretched exponential behavior. However, contrarily to the naive creep 
relation (12) we find that not only U c , but also fj,, depend on temperature. 
Since the phenomenological theory assumes that \i can be computed directly 
from the roughness exponent £ it is important to study the geometrical prop- 
erties of the driven string. For this reason we introduce the averaged structure 
factor, 



S(q) 



dru(r, t)e~ tqr 



(18) 
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The dimensional analysis of this double integral allow us to compute £ from 
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Fig. 10. Structure factor S(q) of the elastic line for the statics. In the presence 
of disorder the long distance behavior is characterized by a roughness exponent 
( eq = 2/3, while in the absence of disorder (see inset) thermal fluctuation gives 
Ct = 1/2. 



S(q) ~ <?~( 1+2 '-'', valid for small q. In Fig. 10 we show the structure factor of 
an elastic string thermally equilibrated at F = 0. We can observe a crossover 
between a short distance regime where thermal fluctuations are dominant 
(C ~ Ct = 1/2) and a long distance disorder dominated regime where we find 
the well known roughness exponent ( ~ £ ei? = 2/3 [55]. For the dynamics 
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Fig. 11. Sketch of S(q) expected for a driven elastic line. 



Dynamics of disordered elastic systems 



15 



(see Fig. 11), when f / one can predict that the short distance behavior 
L < Lx of the clastic string is still dominated by thermal fluctuations (Ct — 
1/2). Note that this thermally dominated regime has nothing to do with the 
regime derived in the previous section and valid up to the scale L opt . In this 
regime disorder is already dominant and barriers are overcomed by thermal 
activation. On the other hand, as already discussed, the finite velocity makes 
the quenched disorder to act as a thermal noise at the largest length scale 
I > L v . Thus, in this case, the expected exponent is also ( v — 1/2 [27]. Finally, 
at intermediate length scales, the physics is determined by the competition 
between disorder and elasticity and characterized by a non trivial random 
manifold roughness exponent. 

A systematic analysis of the v-F characteristics and S(q) show essentially 
two different regimes of creep motion. In Fig. 12(a) and (b) we show the 
structure factor for the two cases. As predicted, we get ( ~ £t = 1/2 for 




-4 =3 -2 =1 -4 =3 -2 =1 

l°g( t l = 27Tn/L) l°g(l = 2TTn/L) 



Fig. 12. After [39]. Structure factor S(q) of the elastic line in the driven case. 
Two regimes are distinguished, (a) In the low disorder/high temperature regime 
the roughness exponent is consistent with ( eq = 2/3. (b) In the high disorder/low 
temperature regime the roughness exponent is clearly bigger than ( eq . 

large q. At a certain scale we observe a crossover between the thermal and the 
random manifold scaling. The location of this crossover decreases as temper- 
ature (disorder) is increased (decreased). We can also observe that the second 
velocity-controlled crossover is not achieved in our finite-size simulation due to 
the very slow dynamics. Interestingly, for the small disorder case (Fig. 12(a)) 
the random manifold scaling gives £ = 0.67 ± 0.05, in excellent agreement 
with the equilibrium value ( eq = 2/3, while a much higher roughness expo- 
nent £ = 0.9 ± 0.05 is found for the strong disorder case (Fig. 12(b)). The 
analysis of the v-F characteristics brings us to the same conclusion: the value 
of the exponent /i is close to the equilibrium value /i = 1/4 for low disorder 
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(high temperature) and departs from this value when the disorder (tempera- 
ture) increases (decreases). In Fig. 13 we summarize all the results. We notice 
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Fig. 13. After [39]. Roughness exponent, and creep exponent, jJ,(T), vs T. The 

dashed line gives the equilibrium roughness exponent £ e<3 — 2/3, and the dotted line 
the purely thermal roughness £t = 1/2. The expected value for the creep exponent 
fj, = 1/4 is also indicated. 

that although the values of £ and /i depart from the equilibrium values, the 
relation (12) seems still to hold, within the error bars for the two exponents. 
This is highly non-trivial since equation (12) is derived from a calculation of 
the barriers in an equilibrium situation. 

We thus have found two regimes of creep motion. The first one occurs when 
the temperature is larger than the strength of the disorder, giving fj, <~ 1/4 and 
( ~ 2/3 as predicted by assuming a quasi-equilibrium nucleation picture of 
the creep motion. This implies that the domain wall has time to re-equilibrate 
between hops, being the underlying assumption behind (12) essentially sat- 
isfied. The second regime occurs for temperatures smaller than the strength 
of the disorder, and is characterized by anomalously large values of both ex- 
ponents. This clearly shows that in this regime the domain wall stays out of 
equilibrium, and that the naive creep hypothesis does not apply. Note that the 
measured roughness exponent is intermediate between the equilibrium value 
and the depinning value Q ep = 1.26 ± 0.01 [24]. The fact that the thermal 
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nucleation which is the limiting process in the creep velocity, is in fact fol- 
lowed by dcpinning like avalanches was noted in the FRG study of the creep 
[29]. Whether such avalanches and the time it would take them to relax to 
equilibrium is at the root of the observed increase of the exponent, is clearly 
an interesting but quite complicated open question. 

6 Conclusions and open questions 

In these short notes we have presented a brief review of the dynamical prop- 
erties of interfaces in a disordered environment. We have in particular focused 
on the response of such interfaces to a very small external force, and the corre- 
sponding very slow motion it entails (so called creep) . Clearly many important 
questions remain to be understood for this problem, for large dimensions the 
microscopic derivation clearly supports the phenomenological one as far as 
the velocity is concerned, but also shows that different length scales enter to 
describe the dynamics. In particular, it predicts a much larger avalanche size 
than initially anticipated. For one dimensional walls the situation is even more 
complex, and the very hypothesis that the wall is constantly in equilibrium 
between two creep processes seems incorrect at least when the disorder is not 
weak enough or if the temperature becomes too low. The deviations such ef- 
fect might entail on the creep exponent is of course important in connection 
with the experimental work. 

Of course these questions are only the tip of the iceberg and more subtle 
questions such as how such a domain wall can age in the presence of the 
disorder are still largely not understood, and more analytical, numerical or 
experimental work is clearly needed to address these issues. 
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